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ABSTRACT 

We present the Allen Telescope Array Twenty-centimeter Survey (ATATS), 
a multi-epoch (12 visits), 690 square degree radio image and catalog at 1.4 GHz. 
The survey is designed to detect rare, very bright transients as well as to verify the 
capabilities of the ATA to form large mosaics. The combined image using data 
from all 12 ATATS epochs has RMS noise a = 3.94 mJy beam~^ and dynamic 
range 180, with a circular beam of 150" FWHM. It contains 4408 sources to a 
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limiting sensitivity of Scr = 20mJy beam~^. We compare the catalog generated 
from this 12-epoch combined image to the NRAO VLA Sky Survey (NVSS), 
a legacy survey at the same frequency, and find that we can measure source 
positions to better than ~ 20". For sources above the ATATS completeness 
limit, the median flux density is 97% of the median value for matched NVSS 
sources, indicative of an accurate overall flux calibration. We examine the effects 
of source confusion due to the effects of differing resolution between ATATS and 
NVSS on our ability to compare flux densities. We detect no transients at flux 
densities greater than 40mJy in comparison with NVSS, and place a 2a upper 
limit on the transient rate for such so urces of 0.004 deg~ ^ . The se results suggest 



that the ^ 1 Jy transients reported by iMatsumura et al.l ( 120091 ) may not be true 



transients, but rather variable sources at their flux density threshold. 
Subject headings: catalogs — radio continuum: galaxies — surveys 



1. Introduction 



The Allen Telescope Array (ATA), a joint project of the Radio Astronomy Laboratory 
of the University of California, Berkeley, and the SETI Institute in Mountain V iew, CA, is a 



radio telescope built on a novel "LNSD" (Large Number of Small Dishes) design ( 1 Welch et al. 
20091 ). It is designed to have a wide field of view (~5 square degrees at 1.4 GHz), enabling 



rapid surveys of large areas of sky. This makes it an ideal instrument to search for transient 
and time-varying radio sources. Since survey speed (the area of sky covered to a given flux 
density limit in a fixed amount of time) is proportional to ND, the product of the number of 
dishes, A^, and their diameter, D ( ISargent fc Welchlll993[ ). this gives the ATA a survey speed 
~ 40% that of the Very Large Array. Plans call for the ATA's expansion from the existing 
42 6m-diameter antennas to a total of 350 antennas, which will result in survey speed an 
order of magnitude faster still than the current configuration. 

The technical design (including the nov el antennas, f eeds, and correlator), and some 
of the early sci ence results, are de s cribed by I Welch et al.l (|2009[ ) . Additional early science 
is described by IWilliams &: Bowerl (j2010[ ) as well as other papers in preparation. In this 
paper we present the ATA Twenty-centimeter Survey (ATATS), which consists of observa- 
tions of 690 square d egrees of sky centered on the ~ 10 square degree Bootes deep field 
( iJannuzi fc Deylll999l ). This is a region with a wealth of multi-wavelength data, including 
data on mid-infrared variability ( iKozlowski et al.ll2010f ). which we will compare with radio 
variability from ATATS in a future paper. 
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ATATS observations were spread over 12 epochs from 2009 January - April (Tabled]). 
Each epoch we aimed to observe as many of the 361 ATATS pointings (§ [3]) as possible; 
in practice, at a given epoch, good data were obtained for between 252 and 318 pointings 
(Tabled]). Observing efficiency and telescope uptime have improved since the ATATS data 
were taken. Occasionally a given pointing failed for some reason - for example, strong 
radio frequency interference (RFI) - and so not all epochs have good data for all pointings 
observed. A few observation attempts resulted in corrupted data or no data at all, due 
to software, hardware, or calibration problems. Epochs with a significant number of failed 
pointings were excluded from the dataset presented in this paper. 

This paper examines the image and catalog made with data from all 12 epochs, and 
compares it to the NRAO VLA Sky Survey (NVSS: ICondon et al.lll998l ). The analysis of the 
images of individual epochs, including transient statistics obtai ned from co r aparing catalogs 
made from data taken at a single epoch, will be the subject of ICroft et al.l (l2010bl ). 



Throughout this paper, we use J2000 coordinates. 



2. Transients and Variables at Radio Wavelengths 

ATATS was designed to verify the performance of the ATA during commissioning (by 
comparison to existing sky survey data) as well as to search for time-varying and transient 
sources. The survey was designed to cover large areas at a comparatively low sensitivity. 

Here we distinguish between sources which are normally above the flux density limit 
(which we term "variable" and "non- variable" , depending on whether or not changes in their 
flux density are detectable from epoch to epoch - although most radio sources are intrinsically 
variable to some degree) and those which have no quiescent radio counterpart (which we term 
"transient"). There is obviously some overlap between transient and variable sources, in that 
some sources which exhibit large outbursts can be seen in quiescent radio emission if observed 
to faint enough flux densities. But it is useful (particularly when comparing catalogs taken 
at different epochs) to distinguish between variable sources which are present in most or all 
epochs (often rotating sources, or sources with variable accretion rates or magnetic fields), 
and transients which appear above the detection limit at one or a few epochs, and are 
undetected at others (often explosive events). Variable and transient sources manifest on a 
range of different timescales in the radio. 



Transi ents have been observe d with timescales from milliseconds (ILorimer et al.l 120071 ) 
to mon ths (IBrunthaler et al.ll2009l ). from a ran ge of progenitors. Thes e include gamma ray 
bursts (IFrail et al.l 120001 ) and radio supernovae (IBrunthaler et al.ll2009l ) - where the progen- 
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itor is destroyed - and flares from a variety of objects, wliicli may or may not also show up 
as variable sou rces during their m ore quiescent periods. Flares are seen from active galactic 
nuclei (AGNs; iFalcke et al.l 119991) at one end of the size sca le, to M- dwarfs (Jackson et al. 



19891 ) , brown dwarfs (iBerger et al.ll2001 



and rotating radio transients (RRATs; 



Hallinan et al. 2007i). p ulsars ( Cognard et al. 1996 ) 



McLaughlin et al.l 120061 ) at the other. 



Some radio transients ( iBower et al.ll2007l : iMatsumura et al.ll2009l ) show no counterparts 



in ver y deep observatio ns at other wavelengths; these may be flares from isolated old neutron 



ical in origin at all - the 



Lorimer et al. 


J2OO7) 


(Burke-Spolaor et al. 



20101), although astrophysical interpretations 



are not yet completely ruled out. 

Variable sources at radio wavelengths are also associated with several diffe rent classes o f 
object and ch aracteristic timescale. Stars often exhib it variable radio emission ( lGiidelll2002l ). 
as do pulsars (lArchibald et al.ll2009l ). microquasars ( jMiller- Jones et al.ll2009l ). and AGNs. 



AG Ns make up alrnost al l of the sources seen in imaging surveys above a few mJy at 
1.4 GHz ( ISeymour et al.ll2008l ). and vary over a range of timescales. Some AGNs exhibit 
IntraDay Variability (IDV) - scintillation of extremely compact components due to ex trinsic 
sources, likely located in our own Galaxy ( iLazio et al.ll2008l : iKoopmans et al.ll2003l ). This 
can cau se flux changes rang ing from a few percent, up to 50% in so-called Extreme Scattering 
Events Jpiedler et al.lll987l ). 



On timescales of months, radio synchrot ron emissio n is sometimes seen t o vary in cor- 
relation with optical (lArshakian et al.l I2OIOI ) or X-ray (IMarscher et al.l 120021 ) emission, as 
uneven accretion rates from the in-spiraling disk propagate through to knots in the outflow- 
ing jets; on timescales of years, even more substantial variations in accretion rate, heating 
of material, and reprocessing of energy released in the accretion disk can occur. 



L ong-term variabil ity of very bright AGNs {S > 1 Jy) has been systematically explored 



[e. g. 



AUer et al.lll992l ). but there has not been systematic characterization of the fainter 



AGN population. Some faint AGNs such as III Zw 2 show extreme variations (IFalcke et al. 
999h. And over pe riods of ~ 10 yr, quasar activity is thought to switch on and off 



(IHopkins et al.ll2005l ). although the duty cycl e is not well unde rstood and accretion may 



stop and start several times over such periods (ISaikia et al.ll2006l ). 



The cadence, frequency, and sensitivity of the observations influences what kind of 
variable and transient objects will be seen in a given survey. We selected an extragalactic 
field, which we observed at 1.4 GHz to relatively low sensitivity, with a cadence of a few days 
to weeks, so almost all variable and non- variable sources seen are AGNs, and any transients 
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would be likely to be radio supernovae or gamma ray bursts, although some sources may be 
members of other classes of object discussed above, or even more exotic objects. 

3. Survey Design and Data Acquisition 

The ATA has 42 dual-polarization antennas, but the ATATS survey was undertaken 
during a period when the array was still being commissioned, and at times, various feeds 
were disconnected, being retrofitted (to improve electrical connections in the feed tip circuit 
board), or otherwise not in operation, so not all antennas, baselines, or polarizations were 
present in all datasets. In any case, at the time of these observations, the correlator was only 
capable of calculating correlations of 64 inputs (either 32 dual-polarization feeds, or some 
other combination of 64 single-polarization inputs from the 42 antennas). Out of a possible 
861 baselines (each with two polarizations) the correlator was therefore capable of producing 
correlations for 496. Due to the fact that not all 64 inputs were populated with good feeds at 
the time of the observations, the mean number of baselines present in the ATATS dataset is 
340, representing an efficiency of 39%. Since the ATATS data were taken, improvements in 
the feeds, and the installation of a second 64-input correlator, have increased the efficiency 
substantially. 

Some of the ATATS data were later flagged due to radio frequency interference (RFI) or 
other problems (§ 14. ip . reducing the observing efficiency still further. Nevertheless, we were 
still able to achieve good throughput in terms of the product of the area surveyed and the 
sensitivity reached, and as noted, the observing efficiency of the ATA continues to improve, 
and has the potential to improve even more dramatically if more antennas are added to the 
array. 

The layout of the ATATS pointings, and the sequence in which they were observed, is 
shown in Fig. [H Pointings were arranged to approximate a hexagonal close-packed grid, 
with the spacing between pointing centers determined by Nyquist sampling, and accounting 
for the convergence of lines of constant right ascension towards the poles. 

Each pointing was imaged once per epoch as a snapshot observation, with an integration 
time of 60s (6 X 10s scans). Every 50 minutes, the calibrator 3C286 was observed for 60s, 
as well as once at the completion of each epoch. Observing always started with pointing 
number 1, and continued until the end of the scheduled observing time, usually around 12 
hours. 

The correlator provided 1024 frequency channels, each of width 102 kHz, across a 
104.9 MHz-wide band centered at 1430 MHz. The 100 channels at each edge of the band 
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were flagged a priori because of poor sensitivity due to filter roll-off, so the effective band- 
width was 84.4 MHz. 



Data Reduction 



As the number of antennas and the volumes of data from radio telescopes increase, 
calibration, imaging, and flagging RFI by hand start to become impractical. As a result, 
automated routines are required to get the best results from the data. Data reduction was 
undertaken using a software suite developed for the reduction o f ATA data, RAPID (Rapid 
Automated Processing and Imaging of Data; iKeating et al.ll2009h . RAPID consists of C shell 
scripts that call MIRIAD ( ISault et al.l 119951 ) tasks. It performs identification and excision 
of RFI, calibrates the visibilities using a primary flux calibrator, and images the individual 
pointings. It then performs phase self-calibration and amplitude flagging. RAPID currently 
operates only on ATA data but is freely available by contacting the authors. Mosaicking, 
catalog creation, and matching to NVSS are performed by additional custom-written scripts. 



4.1. Flagging, Calibration, and Imaging 

RAPID is able to detect and flag RFI on a channel-by-channel basis by looking for 
channels that often rise above the mean amplitude for all cross-correlations. RAPID also 
identifies any problem antennas that show signs of corrupted spectra, and flags them. It then 
makes an image of the calibrator, using full multifrequency synthesis (to avoid bandwidth 
smearing), and iteratively tests for phase and amplitude closure, flags data outside of a 
certain amplitude range, and makes a new image. 

Once all baselines fall within phase RMS and closure limits, iteration stops and an image 
of the calibrator is produced. Baselines that were determined to be bad in the calibrator 
data are also flagged in the source data. RAPID then makes images of the source fields, 
and iteratively applies phase self- calibration and amplitude flagging while attempting to 
maximize the image dynamic range. Once a maximum in dynamic range is achieved, iteration 
stops - for our images, this results in a median number of CLEAN components per image of 
5710, a median RMS of 3.56 mJy beam~^ (2.2 times the theoretical RMS of 1.6 mJy beam~^), 
and a median dynamic range of 120. Images were made using natural weighting which gives 
about a factor 3 lower noise than uniform weighting for the ATA. Superuniform weighting 
was also tried, but the results were unsatisfactory, both in terms of image quality, and much 
higher noise than natural weighting. 



-7- 



o 
in 



o 

o ^ 



o 

CO 



361 342 323 304 285 266 247 228 209 190 171 152 133 114 95 76 57 36 19 
360 341 322 303 284 265 246 227 208 189 170 151 132 113 94 75 56 37 18 



359 340 331 302 283 264 245 226 207 188 169 150 131 112 93 74 55 36 17 



358 339 320 301 262 263 244 225 206 187 168 



130 111 92 73 54 35 16 
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70 51 32 13 



354 335 316 297 278 259 240 221 202 183 164 145 126 107 88 69 50 31 12 



353 334 315 296 277 258 239 220 201 182 163 144 125 106 87 68 49 30 11 
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47 28 9 



350 331 312 293 274 255 236 217 198 179 160 141 122 103 84 65 46 27 8 



349 330 311 292 273 254 235 216 197 178 159 140 121 102 83 64 45 26 7 
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347 328 309 290 271 252 233 214 195 176 167 138 119 100 81 62 43 24 5 



346 327 308 289 270 251 232 213 194 175 156 137 118 



aO 61 42 23 4 



345 326 307 288 269 250 231 312 193 174 155 136 117 98 79 60 41 22 3 
344 325 306 287 268 249 230 211 192 173 154 135 116 97 78 59 40 21 2 

343 324 305 286 267 248 229 210 191 172 153 134 115 96 77 58 39 20 1 
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1.— The field centers of the 361 ATATS pointings. 
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The data were cali brated using 3C 28 6. The cahbrator flux density was determined to 



be 14.679 Jy, using the iBaars et al.l ( 119771 ) coefficients. We did not perform a polarization 
calibration. The ATA produces cross-polarization data, but at the time of this work, the 
calibration and analysis of the cross-polarization data have not been fully tested. 3C 286 is 
~ 10% polarized, but since both X and Y polarizations are present in the data, the effect 
of source polarization on Stokes I flux densities is largely cancelled; ignoring polarization 
calibration will lead to an additional ~ 1% uncertainty in flux density. This is negligible 
compared to other sources of error. 

The flagged visibility data, the bandpasses, and the final images were all inspected by 
eye, but usually RAPID does a sufficiently good job (given optimal input parameters) that 
little human interaction is required. Where artifacts were visible in the final images, or RFI 
or other bad data remained in the visibilities, additional manual flagging was performed. 

For the 12 good datasets an average of 40.6% of visibilities were flagged - this includes 
RFI, the 200 edge channels, and other bad data identified during data reduction. As dis- 
cussed in §131 the ATA was operating at only 39% of its maximum capacity during ATATS, 
so the fraction of data obtained compared to the technical specifications was 23%. This frac- 
tion is now increasing as additional retrofitted feeds are brought online and other technical 
improvements are made, which will improve the fidelity of future surveys such as PiGSS (the 
Pi GHz Sky Survey; Bower et al. in prep.). 

The size of the images made for each pointing was 512 x 512 pixels, where each pixel 
was 35" X 35". The data for all 12 epochs for each pointing were included in each image, 
since the much better {u, v) coverage (Fig. [2]) provides for much better image fidelity. Not 
all pointings had good data at all epochs (for example, because of a failed observation, or 
severe corruption by RFI). Some had bad data at all epochs because they were close to very 
bright sources - in particular, in the regions surrounding 3C286 and 3C295. We included 
only pointings for which at least 9 out of 12 epochs were good - a total of 243 pointings. The 
majority of these (218) had good data for all 12 epochs. The coverage of the final mosaic is 
shown in Fig. [3l 

The images were all made with the same sky coordinate for the reference pixel, to ensure 
that the resulting mosaic should have the correct projection geometry. A typical synthesized 
beam is shown in Fig. |H Since the synthesized beam changes in shape and position angle 
across the mosaic, it was necessary to use an "average" restoring beam to avoid superposing 
different beams from neighboring pointings, which creates problems at the catalog creation 
stage. We examined the beam major and minor axis width as a function of hour angle 
(Fig. [5]), and took the geometric mean of the major and minor axis width at HA = 0, 150", 
as the FWHM of a circular Gaussian which we used as a restoring beam. 
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Table 1. ATATS epochs 



Epoch 


UT Date 


Pointings 


ATAl 


2009 Jan 12 


255 


ATA2 


2009 Jan 19 


271 


ATA3 


2009 Jan 26 


254 


ATA4 


2009 Jan 31 


285 


ATA5 


2009 Feb 7 


285 


ATA6 


2009 Feb 14 


318 


ATA7 


2009 Feb 15 


310 


ATA8 


2009 Mar 27 


261 


ATA9 


2009 Mar 28 


252 


ATAIO 


2009 Mar 29 


270 


ATAll 


2009 Mar 31 


257 


ATA12 


2009 Apr 3 


253 



atat-0072-XX 
atat-007S-YY 



atat-0072-XX 
atat-0072-YY 



u (k\) 



I (kA) 



Fig. 2. — Typical {u,v) coverage for a single pointing at a single epoch (left), and for the 
same pointing when combining the data from all epochs (right). 
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Fig. 3. — The deep field image, made from a combination of all 12 epochs. The solid line 
shows the edge of the mosaic - note that some small regions with fewer than 9 epochs of 
good data were not included in the mosaic. The areas near 3C 295 (the large missing region 
at top center) and 3C286 (the notch cut into the lower right edge) are notable examples, 
but there are also a few smaller regions missing, where > 4 epochs were unusable due to 
RFI or some other problem (see § 14. ip . The greyscale runs from 11.82 mJy beam^^ {3a) to 
118.2 mJy beam-i(30cr). 



Fig. 4. — A typical synthesized beam for ATATS, using data for all 12 epochs for a typical 
pointing (Field 192). The portion of the beam shown here is 75' x 75'. The sidelobes are 
< 10% of the intensity of the main lobe. Each field was CLEANed with its correspond- 
ing synthesized beam, but a common Gaussian restoring beam was used in producing the 
CLEANed images. 



4.2. Primary Beam Correction 



The ATA primary beam h as been measured by variou s methods, including by off-axis 
measurement of poin t sources (MacMahon fc Wrightl l2009l . a method known as "hex- 7";); 



by radio holography (IHarp et al 



2010l ): and by comparison of t he flixx densities of sources 



visible in two or more ATATS mosaic pointings (iHuU et al.ll2010[ ). These measurements give 
consistent results for the FWHM at 1.4 GHz of 150'. We approximated the primary beam 
by a Gaussian model. We also experimented with varying the primary beam FWHM value 
used for mosaicking over a range of values, and found that the canonical value of 150' gave 
us the most accurate and precise fluxes for NVSS sources in our data, and so this value was 
adopted for our analysis. We applied the primary beam correction to the images of each 
good pointing. 

Where the primary beam is well modeled, and where data volumes are small, joint de- 
convolution (where the entire mosaic is imaged and CLEANed at once) is generally desirable 
to maximize image fidelity. ATATS uses snapshot observations (where {u, v) coverage is not 
as good), and the primary beam model used is not as sophisticated, so better quality images 
can be obtained by imaging of individual fields, and subsequent mosaicking. In any case, the 
large amount of ATATS data would make joint deconvolution of the entire mosaic difficult 
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Fig. 5. — The major axis width (squares) and minor axis width (crosses) of the ATATS 
synthesized beam as a function of hour angle, for the observations of 2009 March 31. 
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without some averaging of the input visibihties (which risks degradation of image quahty). 



4.3. Mosaicking 

We made a hnear mosaic of the 243 individual primary-beam corrected images, masking 
each outside a radius corresponding to the nominal half-power points (75'), in order to avoid 
exaggerating the effects of the deviation of the true primary beam fr om a Gaussian at large 



distances from the pointing center (in particular, the 10% sidelobes; iHarp et al.ll2010[ ). 



The final image was transformed to a GLS (global sinusoid) geometry, which has the 
property that each pixel subtends an equal area on the sky. This is important for an image 
covering such a large area, since the curvature of the sky introduces distortions. In projections 
that do not preserve area, the size of the beam would change across the mosaic, resulting in 
errors in the integrated flux densities reported by the source fitting procedure. If we were 
not to perform this transformation on the ATATS mosaic, we would bias the flux densities 
of the sources low by around ~ 2%, and increase the scatter in the measured flux densities 
by ~ 12% when compared to NVSS. 

The RMS noise in this final 12-epoch image, which we term the "master mosaic", is 
3.94 mJy beam~^. The dynamic range, determined by dividing the peak flux density of the 
brighest source in our map (4.14 Jy) by the RMS noise in a 100-pixel box surrounding this 
source (23.0 mJy) is 180. The master mosaic covers 689.63 square degrees of sky - hence, 
the mean master mosaic pixel has contributions from 1.73 pixels of the 243 input 12-epoch 
images, each of which covers 4.91 square degrees. The master mosaic is shown in Fig. [31 and 
a smaller region, which better shows the individual ATATS sources, is shown in Fig. El 



4.4. Catalog creation 

The MIRIAD task SFIND was run on the master mosaic. The RMS was computed over 
a box of side 100 pixels, and sources above 5 times the local RMS were kept. The resulting 
SFIND catalog contained 4843 sources. Of these, 197 were poorly fit by SFIND (indicated by 
asterisks in the SFIND output) and were rejected from the catalog, leaving 4646 sources. Of 
these, 238 sources were in regions where 10% or more of the pixels in a 1° x 1° box centered 
on the source were masked (§ 15. ip . These sources were also rejected, leaving 4408 sources in 
the final catalog (Table [21). 
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RA (J2000) 

Fig. 6. — A 25 square degree region of the deep field shown in Fig. |3l shown at a larger 
scale and slightly different stretch so that the structure of individual sources may be seen. 
All sources with NVSS flux densities brighter than 20mJy (corresponding to 5 times the 
RMS of the ATATS image) are plotted as circles; the size of the circle is proportional to the 
NVSS flux density. The greyscale runs from zero to 118.2 mJy beam~^(30(T). The faintest 
NVSS sources are below the ATATS completeness limit. Hence the fraction of NVSS sources 
without a counterpart in the ATATS catalog increases with decreasing circle size in this 
image. 
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Table 2. The ATATS catalog 



RA n2onoi 


Dec n2000l 


ATATS K^'-'-^yj 






^min ^, ) 


PA (°) 


^^NVSS 


^Nvss y^^^'^y } 


12 55 50.4 


53 03 51 


85.2 


4.9 


164.4 


131.1 


41.8 


1 


79.3 


12 55 52.3 


52 49 59 


102.4 


4.1 


177.4 


137.1 


45.0 


1 


98.1 


12 57 11.2 


53 35 11 


25.1 


5.1 


164.6 


115.7 


28.3 


1 


39.4 


12 57 11.6 


52 43 24 


37.8 


4.2 


204.7 


138.4 


7.7 


1 


34.2 


12 57 37.4 


52 38 26 


105.8 


5.6 


175.1 


136.3 


24.5 


1 


102.4 


12 58 00.1 


52 22 51 


49.8 


5.4 


168.8 


124.8 


11.9 


1 


59.8 


12 58 02.2 


52 51 47 


80.0 


4.2 


194.4 


137.8 


61.9 


1 


70.2 


12 58 14.1 


54 21 53 


678.6 


7.3 


173.9 


142.0 


42.7 


1 


647.0 


12 58 42.5 


54 43 38 


60.6 


4.4 


166.9 


133.8 


44.0 


1 


58.7 


12 59 18.6 


54 32 21 


95.0 


4.9 


190.2 


137.7 


44.9 


1 


82.0 


12 59 37.8 


53 20 44 


40.4 


3.5 


180.2 


167.0 


-29.8 


1 


33.3 


13 00 23.0 


52 35 39 


67.6 


4.7 


170.0 


136.2 


46.1 


1 


71.9 


13 00 38.2 


53 11 49 


37.5 


5.1 


210.2 


155.0 


61.2 


1 


18.1 


13 01 12.0 


53 14 22 


35.6 


5.4 


177.0 


119.3 


-0.9 


1 


40.2 


13 01 22.1 


53 32 56 


120.3 


3.4 


164.2 


141.0 


45.2 


1 


115.5 


13 01 40.5 


54 08 28 


301.3 


3.6 


248.8 


140.9 


65.9 


1 


167.9 
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2 
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1 
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-6.9 


1 


30.2 


13 04 36.5 


54 28 36 


25.1 


4.1 


179.5 
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1 
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1 
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^Number of NVSS matches within 75" of ATATS position 

''Sum of integrated flux densities of NVSS sources within 75" of ATATS position 

Note. — Only a portion of this table is shown here to demonstrate its form and content. A ma<;hine-readable version of the full 
table is available. 
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5. Analysis 

In order to analyze the quality and r eliability of our dat a, we compared the images and 



catalogs to those produced by the NVSS (ICondon et al.lll998l ). a survey at 1.4 GHz (with 45" 



resolution) undertaken with the VLA between 1993 and 1997. One of the goals of ATATS is 
to study source variability; we expect some sources to have changed in brightness, perhaps 
by so much that sources in NVSS may not be present in ATATS or vice versa. However, the 
majority of sources in the NVSS should have similar fluxes and positions when measured 
again by ATATS, and by looking at variations in these quantities we can assess the quality of 
the ATATS data. We can also compare fluxes and positions f rom epoch to epoch in ATATS; 



variability on epoch-to-epoch timescales will be discussed by lCroft et al.l ( l2010bl ). 



We matched the ATATS master mosaic catalog with NVSS in a number of ways, ex- 
perimenting with different source-finding parameters, match radii, and rejection of spurious 
sources in order to obtain the most reliable catalog. We found that a simple circular posi- 
tional match with a match radius of 75" (corresponding to the size of the restoring beam) 
gave a good compromise between too many false matches, and too many true matches that 
might be missed. However, it should be noted that the differing resolutions and sensitivities 
of ATATS and NVSS mean that the "correct" NVSS source is not always matched to the 
corresponding ATATS source. In order to reduce the number of spurious matches of ATATS 
sources to NVSS, we consider only NVSS sources brighter than lOmJy, around 2.5 times 
the RMS of the ATATS master mosaic. We search the ATATS SFIND catalog around the 
positions of NVSS sources, and plot the resulting flux densities of matched sources in Fig. [3 

Sometimes a source unresolved by ATATS is resolved by NVSS into two or more com- 
ponents. This results in apparent changes in flux density for the matched catalog sources, 
since only one of the components (and hence only a part of the object's true flux) at a time 
is matched at these epochs. Because we search around each NVSS source for the nearest 
ATATS source, there are instances where different NVSS sources are each matched to the 
same nearby ATATS source. These cases appear as vertically aligned distributions of points 
in Fig. [71 and are more common at faint fluxes, where the areal density of sources is higher. 
There are also cases where a faint NVSS source matches a bright ATA source because the 
ATA either does not detect the fainter source (and so the NVSS source is matched to a 
nearby brighter, unassociated ATATS source), or because the larger synthesized beam of 
the ATA results in the detection of a single source which NVSS resolves into two or more 
components. In such cases the resulting ATATS flux density is approximately the sum of 
the flux densities of the NVSS components. There may also be cases where the ATA detects 
extended emission which is resolved out by the higher-resolution NVSS. 



Because the flux limit of the NVSS catalog is fainter than that for ATATS, the number 
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Fig. 7. — Comparison of the flux densities of the closest ATATS source (if one exists within 
a radius of 75") to each NVSS catalog position. 
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of NVSS sources within the ATATS field is larger than the number of ATATS sources in the 
same region, and so we can ameliorate the problem of multiple matches of the same ATATS 
source to several NVSS sources by reversing the sense of the match. We searched the ATATS 
catalog for the closest NVSS match, again with a 75" match radius, and plotted the matched 
sources in Fig. [8]). The "tail" of sources to the lower right of Fig. [TJ where a bright ATATS 
source is matched several times to the "wrong" nearby NVSS source, is largely gone. 

However, there are still a number of sources far below the 1:1 line, and fewer sources 
far above. It is equally likely that a source becomes ten times brighter between the NVSS 
and ATATS epochs, as that a source becomes ten times fainter. Were both datasets taken 
with the same instrument, and to the same fiux density limit, the matched catalogs would be 
expected to produce a symmetrical cloud of points around the 1:1 line, with the scatter due 
to a combination of measurement uncertainty and intrinsic variability. The slight asymmetry 
of the point cloud in Fig. [8] is due to the different resolutions and sensitivities of NVSS and 
ATATS. To attempt to account for this, and more reliably compare the two catalogs, we 
"smoothed" the NVSS catalog to match the ATATS resolution, by summing the fluxes of all 
NVSS sources within 75" of the ATATS position. The number of NVSS sources matching 
each ATATS source, and their summed flux density, are shown in Table [21 

The effects of the smoothing can be seen in Fig. [91 The summation of the fluxes of 
NVSS sources within the ATA synthesized beam brings the fluxes more closely into line with 
those from NVSS. However, it is clear that there are still some residual resolution effects, as 
evidenced by the fact that the point cloud is still somewhat asymmetric. 

The spread of the points around the 1:1 line increases towards lower flux densities. This 
is mainly because percentage uncertainties in the flux measurement increase towards lower 
fluxes, but also partly because of intrinsic variability in the sources themselves over the ~ 15 
years between the NVSS and ATATS observations (§ 15. 3p . 

A least squares fit to the data in Fig. [91 has a gradient of 1.006 ± 0.003 and an intercept 
of 1.3 ± 0.5 mJy, so the match of ATATS fiux densities to NVSS is very good. In § 15. 3[ we 
also compare the median fiux densities of ATATS and NVSS and find good agreement. 

5.1. Completeness 

By selecting only ATATS sources which match NVSS, we are able to estimate complete- 
ness of the ATATS catalogs; by selecting those which do not match an NVSS source, we are 
able to identify spurious sources and potential transients. Clearly we wish to reject spurious 
sources while keeping potential transients, and we tweaked the transient detection pipeline 
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Fig. 8. — Comparison of the flux densities of the closest NVSS source (if one exists within a 
radius of 75") to each ATATS catalog position. Note that the sense of the match is reversed 
compared to Fig. [71 
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Fig. 9. — Comparison of the flux densities of sources from the ATATS catalog, with the sum 
of fluxes of aU NVSS sources within 75" of the ATATS positions. 
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to achieve this. We found that the main cause of spurious sources was increased noise at the 
edges of mosaics, so when looking for transient candidates we reject those which are in fields 
with a significant fraction (> 10%) of blanked or masked pixels in a 1° x 1° box centered on 
the source. 

Even without matching the catalogs source by source, we can assess the completeness 
and reliability of ATATS by plotting the flux density distribution of our catalogs and com- 
paring to NVSS. In order to compare source counts over the same area, we culled sources 
from the NVSS catalog in regions where ATATS has no data. The resulting culled NVSS 
catalog, which covers the same area as the ATATS master mosaic, contains 10729 sources 
down to a limiting flux density of 10 mJy. We plot source count histograms of the ATATS 
and culled NVSS catalogs in Fig. [TDl The two histograms are consistent, within the errors, 
until the ATATS counts turn over between 40 and 80 mJy. 

The NVSS histogram truncates abruptly at lOmJy due to the cut we impose on the 
catalog. This is well above the 99% completeness limit of the NVSS, 3.4 mJy, so the NVSS 
counts may be considered a measure of the true source counts for this region of sky. The 
ATATS histogram is marginally higher than that for NVSS for the bins between 320 and 
1280 mJy, and marginally lower for the bins between 40 and 160 mJy (above the points where 
the source counts turn over). This is to be expected since the poorer resolution of ATATS will 
tend to combine the flux of multi-component NVSS sources into a single brighter ATATS 
source, and occasionally to detect extended flux resolved out by NVSS. This should only 
affect a small number of sources though, so the effect on the histogram is not large. 

The ratio of the ATATS to NVSS histograms can also be interpreted as the efficiency 
with which we would detect transient sources of a given brightness, or as a measure of the 
survey completeness. We plot the ratio of the two histograms, which we denote C, in Fig.fTUl 
We can see that ATATS is > 90% complete down to around 40mJy, or approximately 10 
times the RMS noise in the master mosaic, below which the completeness falls off rapidly. 
The shift of some sources into higher flux bins due to the resolution mismatch results in 
some bins with computed C > 1. So some of the sources "missing" at S* > 40mJy are in fact 
detected, but with higher flux density, and the 90% completeness value may be considered a 
lower limit. Matching of individual ATATS sources to NVSS and subsequent examination of 
postage stamp images (§ 15. 4p also supports the conclusion that 90% is a conservative value. 

As shown in Fig. [6l which plots a region of the master mosaic, with NVSS sources 
brighter than 20 mJy overlaid, there is an excellent match between the two surveys, although 
counterparts to some of the faintest sources above the NVSS cut are sometimes marginal or 
non-detections in ATATS. 
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Fig. 10. — Left: Log(N) - log(S) histogram of NVSS sources (black dashed hne) in good 
regions of the ATATS field (i.e. , the culled catalogs - see § 15. ip . and sources detected in 
the ATATS catalog (red solid line) over the same area. Error bars assume Poisson statistics. 
Each bin covers a factor of 2 in integrated flux density, or 0.301 dex. The coarser resolution of 
the ATA tends to combine the flux from a tight cluster of NVSS sources into a single ATATS 
source. Combined with the ATA's additional sensitivity to structures larger than the NVSS 
synthesized beam, this tends to shift some sources from fainter into brighter bins. Right: 
Ratio of the histograms for the ATATS and NVSS catalogs, C = A^atats/^^nyss • This 
represents a measure of the completeness of the ATATS catalogs. Error bars are computed 
by propagating the Poisson errors from the ATATS and NVSS source count histograms. The 
tendency of the ATA to shift sources into higher flux bins results in some bins with C > 1. 
These bins are consistent with 100% completeness (dashed line) within the Poisson errors. 
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We are also able to estimate completeness by adding simulated sources to the ATATS 
mosaic. We took the flux densities of the 10729 sources in the culled NVSS catalog, assigned 
them to random positions within the ATATS mosaic, and created simulated Gaussian sources 
with FWHM 150" at these positions. We then ran SFIND using the same parameters as 
used for the generation of the science catalog (§ I4.4p . The resultant catalog contains both 
the real sources detected in ATATS, and those of the simulated sources that were recovered 
by SFIND. In order to determine our ability to accurately recover the input flux densities, 
we searched for the closest ATATS source (within a maximum radius of 75") to the input 
positions of the simulated NVSS sources. The resulting fiux-flux diagram is shown in Fig.fTTl 
This exercise is similar to the matching of the real NVSS catalog with the ATATS catalog, 
as seen in Fig. [TJ modulo the fact that adding simulated sources increases the number of 
detections in the SFIND catalog from 4843 to 9404 (which will increase source confusion 
somewhat), that resolution effects apply to the match with NVSS, but not to the retrieval of 
simulated sources, and that intrinsic source variability will result in changed fluxes in Fig. [TJ 
whereas no such effect should be seen in Fig. [TTl 

We also searched for ATATS sources within 75" of NVSS sources within the region where 
we have good ATATS data. This provides a similar (and consistent) estimate of completeness 
to the ratio of the ATATS and NVSS histograms in Fig. [TOl 93% of NVSS sources brighter 
than 40mJy are detected by ATATS, as are 97% of sources brighter than 120 mJy. The 
remaining 3% of sources without a match in the ATATS catalog are in fact clearly visible 
in the mosaic images, but are not present in the ATATS catalog due to poor fits or failed 
deb lending (§ 15. 4p . 

5.2. Positional Accuracy 

To determine the positional accuracy of ATATS versus NVSS (which has an intrinsic 
positional accuracy typically < 3") we plot the offsets in right ascension and declination, and 
a histogram of the magnitude of the offset between ATATS and NVSS positions (Fig. [T21) . A 
120" match radius was used for these plots, but as can be seen, the majority of the matches 
fall within a radius corresponding to the ATA synthesized beam (~ 75"), the value adopted 
for the match radius in our analysis. 

The median offset of the ATATS position from the NVSS position is (ARA, ADec) = 
(1"4,0"98), and the RMS scatter in the positional offsets is ((7RA,fTDcc) = (20"4, 14"7). 

As shown in Fig. [6l the positional accuracy is sufficient in most cases to provide a clear 
match between ATATS and NVSS sources, although in some cases a single ATATS source is 
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Fig. 11. — Comparison of the flux densities of simulated sources inserted into the ATATS 
master mosaic, and the fluxes retrieved by SFIND for the nearest source within 75" of the 
input positions. 
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associated with multiple NVSS sources due to the differing resolutions of the two surveys. 



5.3. Variability 

As noted in § 15.11 simulated sources inserted into the image and then recovered in the 
catalog give a measure of the intrinsic scatter due to measurement error and confusion. 
When the ATATS catalog is compared to NVSS, additional scatter is partly due to intrinsic 
variability, which we attempt to quantify here. Fig. [13] shows the same data as in Figs. [7] 
and [TT] (for sources between and 1 Jy), this time on a single plot to aid comparison. Points 
inserted from the simulated catalog have the same resolution as ATATS, and do not cluster, 
unlike when comparing to NVSS, which has higher resolution and multi-component sources. 
As a result, the ATATS-NVSS match has some points (in the lower right quadrant of Fig. [T3l) 
that are due to the ATA detecting a single object that is resolved by NVSS into several fainter 
components. Such sources are not present in the match to the simulated catalog. Even in 
the simulated catalog, however, an inserted source can still be matched erroneously to a 
nearby ATA source; these are the red points in Fig. [13] scattered well below the main cloud 
of points. So even when the resolutions are matched, and sources do not vary with time 
(the red points) there is still quite a lot of scatter just caused by measurement error and 
confusion. Clearly, though, there is some additional spread in the black points, which is at 
least in part due to variability of the sources. 

We use the measured flux densities, S'm^, and comparison flux densites, to compute 
the Pearson product-moment correlation coefficient. 



^ ^ Z]i=l('S'm, Sm){,Sci Sc) 

for the sources shown in Figs. [7] (ATATS near NVSS positions), M (NVSS near ATATS 
positions) and [H] (simulated source retrieval); the values are r = 0.863, 0.969, and 0.973, 
respectively. When we smooth the NVSS catalog (as in Fig. [9]), the correlation is even 
better (r = 0.985), but we cannot apply this smoothing to the simulated catalog because 
there the sources have random positions. A clue to the additional scatter introduced by 
source variability is given by the comparison of r for NVSS near ATATS and for retrieval of 
simulated sources, which shows that measurement uncertainties dominate. 

We can examine this in more detail by computing the scatter in the points for the two 
samples. We define the fractional variability of each data point. 



- 26 - 




Fig. 12. — Left: Positional offsets in RA and Dec for sources detected by botli NVSS and 
ATATS. Right: Histogram of offsets (VARA^ + ADec^) for sources detected by botli NVSS 
and ATATS. Tlie majority of tlie matclies fall within a radius corresponding to the ATA 
synthesized beam (~ 75") - those at larger radii are mostly spurious matches of an NVSS 
source with an unassociated nearby ATA source, where the true match was not detected by 
the ATA. 
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Fig. 13. — The same data as in Fig. [7] (ATATS sources with NVSS matches, plotted in black) 
and Fig. [H] (retrieved flux densities for simulated sources inserted into ATATS, plotted in 
red), shown on a single plot to aid in comparison. 
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and in Fig. HU we plot a histogram in Vi for the ATATS-NVSS match (for the sample as a 
whole, and split into subsamples of sources with Satats < 100 mJy and S'atats > 100 mJy). 
The central peak of the histogram is narrower for the bright sources, showing that their 
fractional variability is lower. All the histograms show the asymmetry towards high Vi due 
to the mismatch of bright ATATS and fainter NVSS sources discussed above. To obtain a 
measure of the true variability we wish to quantify the scatter in Vi without being biased by 
these mismatched outliers, so we choose to measure the interquartile range of V^, as opposed 
to the standard deviation. We compute the interquartile range, IQR^/, in 26 flux density 
bins, with bin edges selected to place an equal number (188) of data points, 5*^,, in each bin. 

By examining the change in IQRy as a function of flux density, we can quantify the 
spread of the points in the fiux-fiux plots. As well as computing IQRy for the ATATS-NVSS 
match, we compute the same statistic comparing the measured flux densities, Sm,, to the 
inserted simulated catalog values, S'c-. We plot IQRy for the two samples in Fig. [T3 Al- 
though there is quite a large amount of scatter in the points, we can see that generally the 
ATATS-NVSS match tends to show a little more variability than the simulated catalog, as 
discussed above, especially above ~ 300 mJy, where the reliability of retrieval of simulated 
sources improves, but the scatter in the match to NVSS is higher. Some resolution effects 
undoubtedly remain, but at least some of this additional variability is presumably intrin- 
sic. The fractional variability generally decreases towards higher flux densities (with the 
exception of the brightest bin, where the mismatched sensitivities and resolutions appear to 
dominate). 

We take the ATATS and smoothed NVSS flux densities (as plotted in Fig. [H]), and plot 
their ratio as a function of ATATS flux density in Fig. [161 Again we see that the scatter 
around the 1:1 line increases towards fainter ATATS flux densities, suggestive of increased 
fractional variability for fainter sources. As discussed in § O an asymmetry in the point 
cloud can be seen, most obviously at flux densities fainter than the ATATS completeness 
limit (to the left of the dotted lines in Fig. [161 Above the completeness limit the point cloud 
is more symmetric about the 1:1 line, although there are still more outliers that appeared to 
get brighter from NVSS to ATATS than those that appeared to get fainter. 

We also plot in Fig. [161 the median flux density ratio in bins of 0.25 dex in S'atats- For 
sources brighter than the completeness limit, the median ratio of S'atats /S'nvss is 0.97±0.10. 
This suggests that ATATS marginally underestimates the fluxes of NVSS sources. This may 
be due to a combination of the differing survey resolutions, the accuracy of the overall flux 
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Fig. 14. — Histogram of the fractional flux variability, Vi (see § 15. 3p for the ATATS-NVSS 
match (black solid line). Also plotted are the histogram for ATATS sources fainter (blue dot- 
dashed line) and brighter (red dotted line) than 100 mJy. Positive values of Vi correspond 
to sources which are brighter in ATATS than in NVSS. 
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Fig. 15. — IQRy, the interquartile range of the fractional flux variability, Vi (see § 15. 3p . 
computed for ATATS sources with NVSS matches (squares) and retrieved flux densities for 
simulated sources inserted into ATATS (crosses). Generally the scatter is larger for the real 
sources. 
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Fig. 16. — Comparison of the flux densities of sources from the ATATS catalog, with the 
sum of flux densities of all NVSS sources within 75" of the ATATS positions. The data are 
the same as those plotted in Fig. [9], but here, in place of the NVSS flux density, we plot the 
ratio of the ATA to NVSS flux densities (black crosses). Sources whose flux density does 
not change from the smoothed NVSS measurement to the ATATS measurement lie on the 
horizontal solid 1:1 line. Sources whose fluxes change by a factor of 2 lie on the horizontal 
dashed 1:2 and 2:1 lines. The vertical dotted line shows the ATATS completeness limit, 
^ATATS = 40mJy. The slanted dotted line shows the same limit for NVSS, Snvss = 40mJy. 
The red solid line shows the median flux ratio in bins of 0.25 dex in Satats, and the associated 
error bars show the interquartile range of the flux density ratio. We plot (as crosses enclosed 
by squares) the 6 sources above the ATATS completeness limit (vertical dotted line) which 
varied by a factor of 2 or more and which appeared relatively compact and isolated in both 
the ATATS and NVSS images (see §[531). 



- 32 - 



density calibration, and CLEAN bias ( ICondon et al.l Il998l ). However, we do not expect 
CLEAN bias to be the dominant effect. The beam sidelobes are < 10% of the main lobe 
(Fig. SJ, the RAPID software is designed not to CLEAN too deeply, and the combination of 
12 snapshot observations into the master mosaic means that the {u, v) coverage is relatively 
good (Fig. ED. 

Some of the sources in Fig. [13] which exhibit more variability than the simulated (non- 
variable) sources are likely real variables, especially those which are significantly brighter in 
NVSS than in ATATS (because the resolution effects discussed above tend to increase scatter 
towards the lower right of the 1:1 line in this plot rather than towards the upper left). Some 
sources that appear highly variable turn out not to be on closer inspection, however, so we 
carefully examined the most extreme cases. 

We examined postage stamps from the ATATS mosaic and from NVSS for the 67 sources 
brighter than the ATATS completeness limit which have measured fluxes a factor of 2 or 
more different (either brighter or fainter) from their flux densities as measured from the 
smoothed NVSS catalog. We judged that only 6 of the sources were relatively compact and 
isolated in both NVSS and ATATS, and may therefore have varied by a factor of 2 or more 
in the 15 years between epochs. These sources are marked by squares in Fig. [161 Postage 
stamps for these sources are shown in Fig. [T71 The remaining 61 candidates had complex 
morphologies or near neighbors (Fig. [T8|) . These are most likely cases where multiple NVSS 
sources were not deblended successfully by the lower-resolution ATATS (see § 15. 4p . although 
we cannot rule out that some may really be highly variable - however, as stated in § [SI we 
would expect such variations to be symmetrically distributed about the 1:1 line if real. 

We plot the candidate highly variable sources on the plot of flux ratio versus flux density 
in Fig. [ini As can be seen, no highly variable sources are seen to get fainter from NVSS 
to ATATS, suggesting that most apparently highly variable sources appear so because of 
resolution effects. Of the 6 good candidates, 2 were brighter than S'nvss = 40 mJy (i. e. , to 
the right of the slanted dotted line in Fig. [T6l) . A search of the NASA Extragalactic Database 
shows that the five brightest sources of the six are known QSOs or radio galaxies (Table [3]). 
We judge that the faintest is also likely a QSO (although not previously classified as such) 
due to its compact radio morphology and variable radio flux. 



5.4. Transient Statistics 



By combining the data from all 12 epochs to make the master mosaic, we should gain 
a factor 3.5 in RMS sensitivity (which scales as the square root of the observing time), but 
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Fig. 17. — The six sources which we judged good candidates (i. e. , those which appear 
compact and isolated, and are brighter than the ATATS completeness limit) to have varied 
by more than a factor 2 in flux density from NVSS to ATATS. Postage stamps are 30' x 30', 
centered on the ATATS position. The 150" circular Gaussian restoring beam is shown in the 
lower right of the ATATS image. The integrated flux density of the ATATS source at the 
center of the image, and the integrated flux density of the closest NVSS source brighter than 
10 mJy within 75", both in mJy, are shown above the ATATS and NVSS images respectively. 
The greyscale runs from —20 to +50mJy beam~^. 




Fig. 18. — Examples of some of the 61 candidates that appeared to have varied by more 
than a factor 2 in flux density from NVSS to ATATS, but have complex morphologies. Such 
sources are prone to resolution effects and less likely to be truly highly variable than are 
isolated sources. Postage stamps are 30' x 30', centered on the ATATS position. The 150" 
circular Gaussian restoring beam is shown in the lower right of the ATATS image. The 
integrated flux density of the ATATS source at the center of the image, and the integrated 
flux density of the closest NVSS source brighter than lOmJy within 75", both in mJy, are 
shown above the ATATS and NVSS images respectively. The greyscale runs from —20 to 
+50 mJy beam~^. 



Table 3. Sources which are hkely to have varied by a factor 2 between NVSS and ATATS 



Catalog ID 


RA (J2000) 


Dec (J2000) 


Satats (mJy) 


CTbg(mJy) 


Snvss (mJy) " 


Probable identification 


ATATS J140138+361117 


14 01 38.8 


36 11 17 


144.8 


5.2 


72.4 


MG2 J140141+3611 (z = 0.507 BL Lac) 


ATATS J144943+451947 


14 49 43.8 


45 19 47 


106.8 


4.4 


49.2 


87GB 144756.5+453159 (z = 1.635 QSO) 


ATATS J134636+290040 


13 46 36.9 


29 00 40 


104.4 


3.2 


36.5 


MG2 J134636+2900 (z = 2.721 QSO) 


ATATS J145242+542824 


14 52 42.6 


54 28 24 


79.7 


4.7 


29.0 


SDSS J145239.17+542808.6 {z = 2.838 QSO) 


ATATS J132817+340339 


13 28 18.0 


34 03 39 


72.5 


4.9 


32.4 


70 1326+3419 


ATATS J133158+284659 


13 31 58.1 


28 46 59 


62.7 


6.0 


23.4 


FIRST J133158.7+28465 



''Sum of integrated flux densities of NVSS sources within 75"of ATATS position 

Note. — The six sources which we judged good candidates (i. e. , those which appear compact and isolated, and are brighter than the ATATS completeness 
limit) to have varied by more than a factor 2 in flux density from NVSS to ATATS. Associations (with redshifts where available) are from the NASA Extragalactic 
Database. 
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in practice we gain somewhat more in image fidelity due to the poor {u, v) coverage of the 
individual snapshots. However, this comes at the expense of sensitivity to short duration 
transients. A transient with a characteristic timescale much shorter than the several days 
between ATATS epochs, and hence which appears only in a single epoch, will have a flux 
density in the master mosaic a factor ~ 12 less than in the image made from the epoch in 
which it occurred. So we expect that by comparing catalogs generated from the 12 single- 
epoch images, we will detect any such transients with ~ 3.5 times better signal-to-noise than 
we are able to by comparing the master mosaic and NVSS. The analysis of the constraints 
on transients and v ariability from epoch to epoch will be the subject of a forthcoming paper 



( ICroft et al.ll2010bf ). For transients with characteristic timescales longer than the 81 days 
between epochs 1 and 12, but shorter than the ~ 15 yr between NVSS and ATATS, the 
master mosaic catalog compared to NVSS (as discussed in this paper) provides the best 
sensitivity. 

To look for such transients, we matched the final ATATS catalog of 4408 sources with 
the culled NVSS catalog (§ 15.11) of 10729 sources, once again with a 75" match radius. 4333 
of the ATATS sources had a match in NVSS. Of the remaining 75 sources without an NVSS 
match, 39 are above 40mJy, the ATATS 90% completeness limit. We examined postage 
stamp images for these 39 sources (from the ATATS mosaic and from NVSS) and concluded 
that none were real transient sources. In the majority of cases, these were sources which 
consisted of multiple NVSS sources (usually as part of a double source - presumably a radio 
galaxy) which were not successfully deblended by ATATS. Examples of some of these sources 
are shown in Fig. [121 

NVSS is complete to a much fainter flux density limit than ATATS, so we can be more 
confident that an ATATS source without an NVSS counterpart is a transient than we could 
be confident that an NVSS source of similar brightness without an ATATS counterpart is a 
transient. Nevertheless, we searched for NVSS sources within the ATATS field without an 
ATATS counterpart. We limited ourselves to the 978 NVSS sources brighter than Snvss = 
120 mJy (i.e., three times the measured ATATS 90% completeness limit) where ATATS 
is over 97% complete (§ 15. ip . We examined postage stamp images of the 26 NVSS sources 
brighter than this limit with no ATA counterpart. All are sources with complex morphologies, 
and in each case it is clear from the images that the NVSS source does in fact have a 
counterpart in ATATS. Although comparing the catalogs using a 75" match radius results 
in only 97% completeness, the fact that fewer than 1 out of 978 NVSS sources brighter than 
120 mJy are truly missed by ATATS when the morphologies of complex sources are examined 
means that we can be > 99.9% confident that no real 120 mJy sources were undetected in 
the ATATS images. 
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We can use the lack of transients seen when comparing ATATS to NVSS to place 
constraints on the population of transient sources brighter than ~ 40 mJy. We see no sources 
in ATATS without a counterpart in NVSS (after visual inspection of the 39 candidates 
without a match within 75"). Poisson statistics give the probability 



Pin) 



nl 



(3) 



of detecting n transients given that r transients are expected over the survey area. 
Poisson 2cr confidence intervals correspond to P{n) = 1 — 0.95 = 0.05, and given that 
we detect n = sources in the field, we derive a 2a upper limit on the transient rate of 
r = —In 0.05 = 3 sources in the 690 square degree survey, or 0.004 deg~^. 

Fig. [20] shows the constraints on transient rates from some of the surveys discussed in §[2l 
compared to the constraint from the comparison o f ATATS and NVSS. I t is hard to reconcile 
the steeply falling source counts suggested by the iGal-Yam et al.l (120061) and ATATS results 



with the detections of 9 transients brighter than 1 Jy reported by iMatsumura et al.l (120091 ) 
unless these transients make up a different population with a very steep cutoff below 1 Jy. 
Th e ATATS 2a upper limit is marginally consistent with the M09 results above 1 Jy, but if 
the iMatsumura et al.l transients have flat or rising source counts towards fluxes < 1 Jy we 
would expect to see them in ATATS. However, it should be noted that the characteristic 
timescale of the transients reported by IMatsumura et al.l is minutes to days, whereas our 
averaging of 12 ATATS epochs means that transients with durations shorter than a few 
days (i. e. , those which might be present in a single ATATS epoch but not in the other 11) 
will have mean flux densities reduced by a factor 12 (and signal-to-noise by a factor 3.5) 
relative to the single-epoch values. Nevertheless, if there are transients which reach average 
flux densities brighter than 12 times our completeness limit (i.e. , ~ 500 mJy and above) for 
at least one minute (the integration time for each individual ATATS epoch), they should 
appear as sources at least as bright as the completeness limit (i. e. , averaged down by up to 
a factor 12 if undetected at other epochs) in the master mosaic. No such sources are seen. 



m 



We will examine t ransient statistics from the images and catalogs from individual epochs 



Croft et all (l2010bf ). 



6. Summary 



Analysis of the combined 12-epoch ATATS dataset shows that the ATA is able to pro- 
duce accurate images of the radio sky, with positions and fluxes matching the legacy NVSS 
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survey at the same frequency. A search for sources seen in ATATS but not present in NVSS 
(within a 75" match radius) turned up 39 candidates above the 40 mJy ATATS completeness 
hmit, but upon inspection, none are plausible transient sources - the majority are double 
sources in NVSS whose components are more than 75" from the ATATS position. This en- 
ables us to place limits on the transient population of fewer than 0.004 transients deg~^, for 
transients brighter than 40 mJy with characteristic timescales > 3 months. The same limit 
applies to transients brighter than ~ 500 mJy with characteristic timescales of minutes or 
longer. Such sources will be discussed further in the next ATATS paper. 

The authors would hke to acknowledge the generous support of the Paul G. Allen Family 
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Fig. 19. — Examples of some of the 39 transient candidates with S'atats > 40mJy and 
'S'nvss = 0. None are truly transients (§ I5.4p . Postage stamps are 30' x 30', centered on the 
ATATS position. The 150" circular Gaussian restoring beam is shown in the lower right of 
the ATATS image. The integrated flux density of the ATATS source at the center of the 
image is shown above the ATATS image. Here there is no NVSS source brighter than 10 mJy 
within 75", so the NVSS flux density is labelled "/ = 0". The greyscale runs from —20 to 
+50mJy beam~^. 
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Fig. 20. — Cumulative t wo-epoch source den sity for radio transients as a function of flux 
density, based on Fig. 9 of iBower et al.l (120071 ). The solid black line shows the rate measured 
for transients with characteristic timescales < 7 days from lBower et al.l (120071 1 . while the bro- 
ken black lines show the 2a upper and lower bounds. The red das hed lines show 5""^'^ and 
S*"^'^ curves. The arrows show 2a upper limits for transients from iBower et al.l with a 1 yr 
timescale (labeled B07.1) and a 2 month timescale (labeled B07.2), and for t ransie nts from 
the comparison of the 1.4 GHz NV SS and FIRST surveys by iGal-Yam et al.l (|2006|). la beled 



Frail et al. 



(120031 ) sur- 



G06; from the ICarilli et al.l (|2003|) survey (labele d C03) ; from the 

vey (labeled F03); and from the iMatsumura et al. ( 2009 ) survey (l abeled M09). The upper 
limits from ATATS (§ 15. 4p are only marginally consistent with the IMatsumura et al.l detec- 



tions. C onstraints from data from individual ATATS epochs will be discussed in 
(l2010bh . 
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